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ABSTRACT 

Small departures from a homogeneous isotropic spacetime create observable features 
in the large-scale anisotropy of the cosmic microwave background. We cross¬ 
correlate the maps of the cosmic microwave background anisotropy from the Cosmic 
Background Explorer (COBE) Differential Microwave Radiometers (DMR) 4-year 
data set with template maps from Bianchi VIR cosmological models to limit global 
rotation or shear in the early universe. On the largest scales, spacetime is well 
described by the Friedmann-Robertson-Walker metric, with departures from isotropy 
about each spatial point limited to shear a/ Hq < 10“® and rotation uj/Hq < 6 x 10“® 
for 0.1<Do = l- 
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I INTRODUCTION 


One of the goals of cosmology is to understand the large scale structure of the universe. 
Gravitational instability models hold that large-scale structure forms as the result of 
gravitational amplihcation of initially small perturbations in the primordial density 
distribution. The initial “seeds” may result from quantum fluctuations in a scalar 
held (inhation) or causal ordering following a phase transition with broken symmetry 
(topological defects), but in all models the evolution occurs in a “background” 
cosmology described by the Friedmann-Robertson-Walker (FRW) metric. 

The anisotropy of the cosmic microwave background (CMB) provides an 
observational test of the the assumption that on large scales spacetime asymptotically 
approaches the FRW metric. The detection of huctuations in the CMB provides 
support for a nearly homogeneous spacetime, limiting density perturbations in the 
early universe to the level 5p/p ~ 10“®. The CMB anisotropy also provides a direct 
test of the assumption of isotropy about each spatial point. Small deviations from 
the FRW metric lead to observable signatures in the CMB. Open or hat models 
with global rotation or shear will exhibit a spiral pattern of temperature anisotropy 
resulting from the handedness of the geodesics propagating through an anisotropic 
spacetime. Further geodesic focusing creates “hot spots” in open models (density 
Co < 1), while closed models exhibit a pure quadrupole pattern. 

Several authors have used the CMB to limit rotation uj and shear a in the 
universe BBBi- Prior to the detection of CMB anisotropy, Collins & Hawking and 
Barrow, Juskiewicz, & Sonoda used upper limits on the CMB quadrupole amplitude 
to limit uj/Hq < 10“® for hat or moderately open (Cq ~ 0.3) models. Smoot used the 
quadrupole detection from the hrst-year COBEDMR sky maps to limit uj/ Hq < 10“®. 
Recently, Bunn, Ferreira, & Silk htted the full spiral pattern from models with global 
rotation to the 4-year DMR data and derived limits uj/Hq < 3 x 10“^. 

The DMR anisotropy data are dominated by a power-law spectrum of 
huctuations which are ill-described by the pattern of angular anisotropy predicted 
for models with global rotation or shear. Limits to rotation or shear based on full- 
sky “template” maps must account for the presence of this power-law component. 
Bunn, Ferreira, & Silk use a least-squares ht including only instrument noise, and 
then account for chance alignment of the spiral template map with features in the 
power-law component using Monte Carlo simulations. In this paper, we present an 
independent analysis using a formalism which expressly includes the dominant power- 
law component in the htting process, and derive upper limits uj/Hq < 6 x 10“^ for 
0.1<Do = l- 

II GEOMETRY 

The FRW geometry is a special case of more general solutions to the Einstein held 
equations. Relaxing the requirement of isotropy about each point leads to more 
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complicated solutions which contain the FRW metric as a special case. The Bianchi 
models Vllh, VIIq, and IX are the most general solutions for a homogeneous but 
anisotropic universe which contain the open, flat, or closed FRW metric as a special 
case. The metric is dehned as 


= -n^n„M + ^(abE^E® 


where n^ = —t-^ is the normal to the spacelike surfaces of homogeneity, is a 3 x 3 
matrix depending only on t, and E^ are the three invariant covector helds in the 
surfaces of homogeneity such that 





e;^e 


c 

U 


where Cgc are the structure constants (see, e.g., Ellis & MacCallum P). The matrix 
qab can be expanded into a part representing the expansion and a trace-free part 
representing the anisotropy, 

5'ab = e^“[e^^]AB- 

The metric anisotropy in these models is characterized by the vorticity uj and the 
shear a (differential expansion in orthogonal directions). 




^ [(exp(^))*(exp(-/?)) + (exp(-^))(exp(/?))*]ij 

where u is the four-velocity, is an anti symmetric tensor and the dot represents 
a time derivative. The FRW metric is recovered in the limit uj/Hq —>■ 0 and a/Ho —>■ 0 
where Hq is the Hubble constant. 

Photons propagate along geodesics from the surface of last scattering to the 
observer; for small metric anisotropies the CMB anisotropy is 

AT ■ 

-jr = (p'u')r - (p'u')e - 


where p‘ are the direction cosines of the null geodesic and the subscripts E and R refer 
to emission and reception, respectively. The hrst two terms are the Doppler anisotropy 
from the motion of the receiver and the surface of last scattering, while the last term 
represents the geodesics integrated over the anisotropic metric. Various authors have 
solved the geodesic equations for the observed pattern of CMB anisotropy ll|, @, H . In 
this paper, we specialize to Bianchi VIR model (C|^ = C 21 = 1 and = C|^ = y/h), 
which contains the Bianchi I, V, and VIIq models in the appropriate limit. Two effects 
are important. The intrinsic handedness of the tensor (Jij causes the geodesics to spiral 
in open or flat models. The resulting temperature pattern is given by 

^ J ^(^R)sin0R + S(6 »r)cos0r] (1) 
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(Barrow, Juskiewicz, and Sonoda 0), where (6'r, 0r) specify the photon direction]^, 


dl(0R) 


Ci[sin 6 *r — C 2 (cos' 0 ^ + Sv^sin-^^;)] 
s(l — s^) sin-^dr 


rm 

+ C 3 


te (1 + s 2 ) 2 sinh"^(v^ r/ 2 ) ’ 


( 2 ) 


-B(6*r) = Ci[3'\/hsin6*R — C'2(sin'0£; + Sv^cos'^r;)] 

rni s(l — s^) cos-^dr 




Jte (1 + s^)^sinh"^(-\/h r/ 2 ) ’ 
s = tan exp[-A/h(r - tr)], 


^ = (t - 'Tr) - hi jsin^ ^ + exp[2A/h(r - tr)] cos^ 

and dt = e“dr. The constants Ci, C 2 , and C 3 are dehned as 

1 


^3 = 


^ 3f2o^ 

2se(1 + ze) 
1 + s|; 
4^(1 - 


C 2 = 




where 2 ; is redshift and the parameter 


(3) 

(4) 

(5) 


( 6 ) 


determines the scale over which the principal axes of shear and rotation change 
orientation. We obtain the Bianchi VIIq model in the limit ff ^ 1, h —> 0 with 
X hnite; we obtain the Bianchi I model as x —cxd and the Bianchi V model as 
h —>• cx). 

The ratio u/Hq determines the amplitude of the anisotropy, while the 
parameter x determines the pitch angle of the spiral. A signihcant amount of the 
resulting CMB anisotropy is present at high-order moments than the quadrupole. A 
second effect is geodesic focusing, present only in open models. The component of 
the geodesics along the symmetry axis transforms in these models as 


tan 




exp(rR 



^ The direction an observer sees on the sky is given by 0 = tt — 0r, 


(f) = 'n + 4>Yi- 
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Since 


sinh 



- 1 
l + z ’ 


the spiral anisotropy discussed above is focused into a single hot “navel” ~ hlo 
radians in scale oriented along the symmetry axis of the metric. Figure shows a 
polar projection of the anisotropy for hlo = 0.3 and x = 0.5. 


Ill ANALYSIS 

The DMR 4-year sky maps are well described by a superposition of instrument noise 
and nearly scale-invariant CMB anisotropy, in the sense that the observed sky falls 
near the median of the distribution of simulated skies which include only these two 
components 0 . We test for the presence of a third component of unknown amplitude 
whose spatial distribution is described by a Bianchi VIR model, 

where is the thermodynamic temperature of the scale-invariant CMB 

anisotropy, AX is a “template” map of a specihc Bianchi model, and n is the DMR 
instrument noise. We do not analyze the closed Bianchi IX model, which exhibits 
neither spiral anisotropy nor geodesic focusing. We use a weighted combination of the 
6 channel maps from the DMR 4-year data, from which a model of Galactic emission 



North Hemisphere 


South Hemisphere 


Figure 1: Anisotropy template map (arbitrary units) for a Bianchi VIR model with 
Do = 0.3 and x = 0.5. Note the differences between the hemispheres from geodesic 
focusing. 
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has been removed (the “correlation” map of Refs, p, |^). To farther rednce residual 
Galactic contamination, we analyze only the high-latitude portion of the resulting 
map, here dehned as latitude |6| > 20° with custom cutouts at Ophiuchus and Orion 
0 . For computational efficiency, we degrade the maps one step in pixel resolution. 


leaving 954 high-latitude pixels, each of size 5?2. 

We estimate the correlation coefficient a by minimizing 


x" = i: (r - aX), 

a,b 


iM-^)ab (T - aX),, 


(7) 


where T is a vector consisting of the DMR temperatures in each high-latitude pixel, 
X is a similar vector for the template map, and M is the covariance matrix between 
the elements of T. Several authors have demonstrated that can adequately 

be described as a Gaussian random field with a scale-free power spectrum [0. We 
assume that the Gaussian component is uncorrelated with the geodesic component 
AX. The covariance matrix (including instrument noise) thus becomes 


Mab = {TaT,) = J2{2i + l)WfCePe{na ■ n,) + 


iV„ 




( 8 ) 


and is a function of the instrument noise and power spectrum C^, where Wj is the 
experimental window function that includes the effects of beam smoothing and finite 
pixel size, Piiua ■ Ub) is the Legendre polynomial of order i, lia is the unit vector 
towards the center of pixel a, and Na is the number of observations for pixel a. In the 
limit that only instrument noise is considered {Qrms-ps = 0), the covariance matrix 
M is diagonal and Eq. (|^) reduces to a noise-weighted least-squares estimate of a. 

We adopt a scale-invariant power spectrum |]Tl| , 

r(£+(n-l)/2)r((9-n)/2) 


C(,{,Qrms—PSi'^^ — l^')Qrms—PS 


r(£+ (5-n)/2)r((3 + n)/2)’ 


(9) 


with amplitude Qrms-PS = 18 /rK and index n = 1, and set Gq = Gi = cxo to 
assign zero weight to the monopole and dipole terms (in practice, we use Co = 
Cl = 10® /iK^). Strictly speaking, the power spectrum in Eq. (^ represents a flat 
cosmology, while the Bianchi VIR template maps represent open models. Although 
our choice of an exact scale-invariant flat model may not represent the true GMB 
distribution in detail, it is certainly adequate for these purposes: a variety of analyses 
have shown that the DMR data can not discriminate at high confidence between a 
scale-invariant power spectrum and alternate (e.g., open) models |^. Our results 


are not dependent on the detailed model normalization: using Qrms-PS = 24 /iK and 
n = 0 changes our limits by only 20%. 

Equation m provides an analytic solution for the coupling coefficient a 
between the template map X and the DMR data. 


a = 


Y.a,b Tg (M )abXb 

J2a,b (M“^)aftXb 


( 10 ) 
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assuming a specific orientation of the two maps. In practice, the relative orientation 
of the template is unknown. For a selected Bianchi model specihed by the parameters 
[flo, a:], we create a template map by numerical integration of Eqs. (0) - (0), and hnd 
the htted value of the correlation coefficient a and its uncertainty 6a in each of 13824 
different orientations on an approximately 10° angular grid: we step the symmetry 
axis to point to each of the 384 unit vectors on a quadrilateralized sky cube 0 
at resolution index 4, then rotate the template by 10° steps in azimuth about each 


new symmetry axis. The uncertainty 6a = J2a,b varies slightly 

with angular orientation as different features in the template map X move into and 
out of the low-latitude region excluded by the Galactic cut. Accordingly, we use the 
parameter F = a/6a to assess the statistical signihcance of the correlation in each 
specihc orientation, and select the largest |F| to denote the “best” orientation. 

Even in the absence of geodesic effects in the CMB, chance alignments between 
the Bianchi template map and the power-law component of CMB anisotropy can 
create non-zero correlations. If the relative orientation of the DMR and template 
maps is held constant, F will follow a normal distribution with zero mean and unit 
variance over an ensemble of CMB maps. By selecting the largest |F| from a large 
number of possible orientations in a single CMB map with a fixed pattern of noise and 
anisotropy, we maximize |F| over an ensemble of chance alignments and thus expect 
the largest |F| to lie on the tail of a statistical distribution. 

We test the accuracy of the htting process and assess the statistical signihcance 
of the results using Monte Carlo simulations. For each choice of Bianchi parameters 
[Co, 3;], we synthesize 100 realizations of the power-law component of CMB anisotropy 



Figure 2: Histogram of parameter F = a/6a with input a = 0. Chance alignments 
produce a “best htted” value between 3 and 5 standard deviations from zero. 
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[Eq. (^], add instrument noise following the DMR 4-year observation pattern, and 
derive the best fitted correlation coefficient a from the maximum of |r| as described 
above. The distribution in fitted orientation is fiat: there is no preferred direction 
for the symmetry axis as might be caused by the Galactic cut. Figure || shows the 
distribution in |r|. Most of the realizations have |r| 3.7; 95% of the realizations 

have |r| < 4.5, independent of the parameters [Go,a:]. 

We also synthesize additional realizations which include a component OqX 
with r = 7 to be safely above the noise. When the input template X is exactly 
aligned on our angular grid, the distribution of htted a has mean (a) = ao, and 
correctly selects the input orientation. When the input template is misaligned, we 
recover a slightly reduced amplitude (a) = 0.9q;o with orientation correctly centered 
on the closest grid point. 


IV RESULTS 


We hnd no statistically signihcant correlations between the DMR 4-year CMB map 
and any Bianchi VIR model with O.lfDM and 0.17x710. The most signihcant 
correlation occurs at Dq = 0.15, x = 0.2, for which T = 4.37. A coefficient this 
large occurs in 7% of the null simulations. The minimum value occurs at Dq = 0.95, 
X = 2.0, for which T = 2.91 (1% of simulations). The range of htted coefficients a 
over the grid of values [Do, 3;] agrees well with the distribution expected from chance 
alignments of instrument noise and power-law CMB anisotropy: we hnd no compelling 
evidence for the presence of a third component traced by any Bianchi VIR model. 
We thus adopt the value T = 4.5 at each grid point in [Do,x] as a 95% conhdence 
upper limit to the amplitude of microwave anisotropy traced by the corresponding 
Bianchi template. 

Figure ^ shows the resulting limits to shear a and rotation u. In the limit 
X ^ oo, the template pattern approaches a pure quadrupole and (t/Hq is limited by 
confusion with the observed quadrupole amplitude and quadrupolar emission from the 
Galaxy. For x < 2 the spiral pattern dominates and we obtain a limit a/Ho < 10“®, 
corresponding to rotation uj / Hq < b x 10“® for hat models and uj/Hq < 6 x 10“® for 
open models. 

Our limits on a/Ho are approximately a factor of three below those of Bunn, 
Ferreira, & Silk. This results primarily from our inclusion of the power-law component 
of the CMB anisotropy in the covariance matrix M, which accounts for the chance 
alignments expected (in an ensemble average sense) between the template and power- 
law components. Bunn, Ferreira, & Silk use a noise-weighted ht. 


T.xf/af 


( 11 ) 


(c/Eq. (7) in Ref. 0), and set limits on cosmological parameters by comparing the 
value from the DMR data to the distribution of values from Monte Carlo simulations 



with a = 0. We may reproduce this technique in our analysis by keeping only the 
noise term in Eq. (|[) and examining the resulting distribution in |r| for Monte Carlo 
simulations that include the power-law contribution in the sky maps. If the power- 
law contributions are neglected in the fitting process, the distribution of |r| shifts 
to larger values and becomes substantially broader, resulting in weaker limits to the 
rotation and shear. If we repeat our results using only the noise term in Eq. (H), we 
obtain weaker limits compatible with those of Bunn, Ferreira, & Silk. 

We conclude that the observed microwave sky shows no evidence for geodesic 
effects caused by an anisotropic spacetime. On the largest scales, spacetime is well 
described by the FRW metric, with deviations from homogeneity 5p/p ~ 10“® and 
departures from isotropy a/H q < 10“®. 
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Figure 3: 95% conhdence level upper limits to shear and rotation from the DMR 
4-year CMB map. (top) Shear a/Ho x 10^°. (bottom) Rotation uj/Hq x 10^°. 
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